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Abstract 

This article discusses spin transport in systems with spin-orbit inter- 
actions and how it can be understood in a semiclassical picture. I will first 
present a semiclassical wave-packet description of spin transport, which 
explains how the microscopic motion of carriers gives rise to a spin current. 
Due to spin non-conservation the definition of the spin current has some 
arbitrariness. In the second part I will briefly review the physics from a 
density matrix point of view, which makes clear the relationship between 
spin transport and spin precession and the important role of scattering. 
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Glossary 

• Extrinsic effect an effect which has an explicit dependence on the form or 
strength of the disorder potential. 
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• Intrinsic effect an effect which does not depend expHcitly on the form and 
strength of the disorder potential. 

• Semiclassical theory a theory in which a particle's position and momentum 
are considered simultaneously. 

• Spin- orbit interaction a relativistic interaction between the spin of a par- 
ticle and its momentum (which is associated with its orbital motion.) 

• Steady-state spin current a flow of spins induced by an electric field. 

• Steady-state spin density a net spin density induced by an electric field. 

1 Definition of the subject and its importance 

Spin transport refers to the physical movement of spins across a sample and, 
if spin were a conserved quantity, one could make a straightforward distinction 
between spin-up and spin-down charge currents. The recent upsurge of inter- 
est in spin transport is, however, motivated by systems in which spin is not 
conserved due to the presence of spin-orbit interactions, which give rise to spin 
precession. Here, due to non-conservation of spin the spin current is not well 
defined [H [H [3] . Spin transport in these cases usually does not involve charge 
transport as the charge currents in the direction of spin flow cancel out. Finally, 
in certain materials, spin currents are accompanied by steady-state spin densi- 
ties. The appearance of a spin density is not a transport phenomenon, but it is 
a steady-state process and is intimately connected to spin transport. 

The word semiclassical as used in this work refers to theories which consider 
the position and momentum of a particle simultaneously. Semiclassical pictures 
are intuitive and useful in descriptions of transport, particularly in inhomoge- 
neous systems and in spatially dependent fields, which typically vary on length 
scales much larger than atomic size. 

In recent years, steady progress has been made towards realization of conve- 
nient semiconducting ferromagnets and spin injection into semiconductors from 
ferromagnetic metals [11 [SJ [SI [7] yet spin injection from a ferromagnetic metal 
into a semiconductor is hampered by the resistivity mismatch between the two 
|5]. This is one factor, in addition to basic science, motivating the search for an 
understanding of the way spins are manipulated electrically. The last few years 
have seen many experimental advances in spin transport, and spin currents have 
been measured directly |9l[10] and indirectly [TT l nt l [T3| . fTl [T3] . 

2 Introduction 

Novel physical phenomena that may lead to improved memory devices and ad- 
vances in quantum information processing are closely related to spin-orbit in- 
teractions. [16j Spin-orbit interactions are present in the band structure and 
in potentials due to impurity distributions. Spin-orbit coupling is in principle 
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always present in impurity potentials and gives rise to skew scattering. Band 
structure spin-orbit coupling may arise from the inversion asymmetry of the 
underlying crystal lattice [17| (bulk inversion asymmetry), from the inversion 
asymmetry of the confining potential in two dimensions [18] (structure inversion 
asymmetry), and may be present also in inversion symmetric systems. |19) 

Although many observations in this entry are general, the discussion will fo- 
cus on non-interacting spin-1/2 electron systems, which are pedagogically easier. 
The Hamiltonian of these systems typically contains a kinetic energy term and 
a spin-orbit coupling term, = + where m* is the electron ef- 
fective mass. In spin-1/2 electron systems, band structure spin-orbit coupling 
can always be represented as a Zeeman-like interaction of the spin with a wave 
vector-dependent effective magnetic field ftk, thus — (1/2) cr flk- Common 
examples of effective fields are the Rashba spin-orbit interaction, jTH] which is 
often dominant in quantum wells with inversion asymmetry, and the Dressel- 
haus spin-orbit interaction, [17] which is due to bulk inversion asymmetry. The 
spin operator is given by s"' = {h/2) a" , where g" is a Pauh spin matrix. The 
spin current operator in these systems will be taken to be J" = (1/2) {s'^,-u*}, 
where the velocity operator is = (l/?i) dHh/dki. 

An electron spin at wave vector k precesses about the effective field fJfc with 
frequency Jlfc/Zi = \Q,k\/h and is scattered to a different wave vector within 
a characteristic momentum scattering time Tp. I will assume in this work that 
EpTplh ^ 1, where is the Fermi energy, which is equivalent to the assumption 
that the carrier mean free path is much larger that the de Broglie wavelength. 
Within this range, the relative magnitude of the spin precession frequency f2fc 
and inverse scattering time I/t^ define three qualitatively different regimes. In 
the ballistic (clean) regime no scattering occurs and the temperature tends to 
absolute zero, so that £F7p — * oo and ^IkTp/h oo. The weak scattering regime 
is characterized by fast spin precession and little momentum scattering due to, 
e.g., a slight increase in temperature, yielding epTp/h ^kTp/h 1. In the 
strong momentum scattering regime epTp/h 3> 1 S> ilfc Tp/h. I will concentrate 
on effects originating in the band structure, the observation of which requires 
the assumption that the materials under study are in the weak momentum 
scattering regime. Electric fields will be assumed uniform. 

The first part of this article will present a semiclassical theory of spin trans- 
port, identifying the terms responsible for spin currents in the microscopic dy- 
namics of carriers. Spin non-conservation as a result of spin precession leads 
to several possible definitions of the spin current, which emerge out of the spin 
equation of continuity. The second part presents a different point of view, which 
explains aspects not easily captured in the semiclassical approach. The steady- 
state density matrix is shown to contain a contribution due to precessing spins 
and one due to conserved spins. Steady state corrections ex Tp are associated 
with the absence of spin precession and give rise to spin densities in external 
fields. [20l[2ll[22l[2l[2l[l6l[25l|22l[28] Steady state corrections independent of 
Tp are associated with spin precession and give rise to spin currents in external 

fields. [iiiiiiaiiniiiiiiiiiiiiiiiiaisiisniiniiMiMiiMinaiMisB 
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[MlliniSIlllUlSlllliSliS] Scattering between these two distributions induces 
significant corrections to steady-state spin currents. 



3 Spin currents in electric fields 

3.1 Wave-packet picture of spin transport 

This section presents a semiclassical theory of spin transport valid for a general 
spin-orbit system. The semiclassical method is a suitable approach to the study 
of transport, because, typically, in the relevant systems the external fields vary 
smoothly on atomic length scales. All information about the system is taken to 
be contained in the band structure, thus allowing a description of spin transport 
which does not make reference to the detailed form of the spin-orbit interaction. 

The system under study is regarded collection of carriers, whose 

semiclassical dynamics in a non-degenerate band i are described by a wave 
packet [17], with its charge centroid having coordinates (rc,kc) 

\w,)^ /d3fca(k,t)e^''>,(r„k,i)). (1) 



In the above, the function a(k, t) is a narrow distribution sharply peaked at k^, 
the phase of which specifies the center of charge position r^, while \ui{rc,'k,t)) 
are lattice-periodic Bloch wave functions. The size of the wave packet in mo- 
mentum space must be considerably smaller than that of the Brillouin zone. In 
real space, this implies that the wave packet must stretch over many unit cells. 

The external electric field drives the center of the wave packet in fc-space 
according to the semiclassical equations of motion 

ds 

nr, = ^ - X n, (2) 
hkc = gE, 

with q the charge of the carriers, Si the band energy, and the Berry curvature 

dui dui 
^^ = *^9k"""5k^- 

The electric field also gives rise to an adiabatic correction to the wave functions, 
which mixes the states making up the wave packet. The wave functions \ui) 
therefore have the following form: 

where the (/<i are the unperturbed Bloch eigenstates. The \u.i) form a complete 
set and retain the Bloch periodicity. 
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Figure 1: For a particle of finite extent the charge and spin distributions in real 
space are in general do not coincide. The same is true of the charge and spin 
distributions in reciprocal space. 



The distribution of carriers is described by a function /. When scattering is 
present, the distribution function satisfies the following equation: 

df . df ■ df fdf\ 

where {^)coii is the usual collision term. In independent bands, in the relax- 
ation time approximation, the collision term takes the form , with /o the 
equilibrium distribution and Tp the momentum relaxation time. In the Boltz- 
mann theory, the change in the distribution function with time arises through 
the drift terms, which are determined from the semiclassical equations of mo- 
tion, as well as through scattering with other carriers, with localized impurities 
or with phonons. For transport in a non-degenerate band, it is consistent to 
ignore interband scattering effects in the weak scattering limit. In this case the 
relaxation time is a scalar quantity. The effects of interband coherence due to 
scattering will be explored in the next section. 

In order to obtain expressions for macroscopic quantities of interest, such as 
densities and currents, one needs to carry out a coarse graining by averaging 
over microscopic fluctuations. In classical dynamics this coarse graining is per- 
formed by means of a sampling function, which is smooth and has a significant 
magnitude only in a finite range [48j . This range is large compared to atomic 
dimensions, but small compared to the scale of variation of the distribution func- 
tion. Moreover, it has a rapidly converging Taylor expansion over distances of 
atomic dimensions, and its form does not need to be specified. This method has 
a close analog in wavepacket dynamics, where the sampling function is replaced 
by a (5- function. 

It is crucial to recognize that, in general, the center of spin and the center of 
charge are distinct, since the wave packet samples a range of wave vectors and 
the spin is usually a function of k. Following the line of thought outlined above. 
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Figure 2: In the presence of spin-orbit interactions the spin distribution of a 
particle changes in time. The horizontal axis may represent position or wave 
vector. 



the spin density is defined to be (henceforth kc will be abbreviated to k) 

S"{R,t) = f [ d^kd\,f{r,,k,t){S(R-r)s^), (6) 



where the bracket indicates quantum mechanical averaging over the wave packet 
with charge centroid (rc,k). As the 5-function has operator arguments, it will 
be regarded as a sampling operator, whose expectation value yields a spatial 
average, evaluated at position r. To account for the fact that spin is not con- 
served, a new quantity is introduced, which will be referred to as the torque 
density, defined by 

T-(R,i)= J J d^kd\j{r,,k,t){6{-R-r)f''). (7) 

T°' in the above stands for the rate of change of the spin operator, given by 
j-[H, 5°^], and symmetrization of products of non-commuting operators has been 
assumed. Finally, the microscopic spin current density is defined as: 

J""(R,t) = j j d'kd\,f{vM){5{R~r)s''v). (8) 

We obtain the following continuity equation for the spin density and current: 

V-J-"=T"+^". (9) 



dt 

The equation of continuity contains a bulk source term, which coincides with 
the torque density and acts as a mechanism for spin generation. Similar source 
terms are associated with nonconserved quantities, for example, in quantum 
electrodynamics and in Maxwell's equations. The last term in ^ represents 
the scattering contribution, which will be discussed further below. 

Let us discuss the terms in the equation of continuity, beginning with the 
spin density. The argument of the sampling operator can be expressed as [r — 
Tc — (f — Tc)], and, as the second term is of atomic dimensions, the sampling 
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operator can be written as a Taylor expansion about (f — Tc). The density can 
therefore be re-exprcsscd, in terms of macroscopic quantities, as 



S''{R,t) = p'^{R,t) - Vr • P'^R,t), (10) 

where summation over repeated indices has been assumed. In the above, the 
monopole density is given by 

p'^'R,t) = J yd3fcdV/(rc,k,t)(r)5(R-re) = J d^kf{s^)\r^=n, (11) 

where / in the second hne, and henceforth, is to be understood as /(R, k, f), 
and the dipole density is 

P^"(R,i)= J J d3A;d3rc/((r-re)s")5(R-re) = J d"" kf p'^lr^^n- (12) 

The average spin of the wave packet has been denoted by (s^), and the spin- 
dipole is defined to be p'^'^ = ((f — rc).s'^) |r,.=R- It will be seen that the first 
term in the density is the average of a monopole density located at Tc, while 
the dipole term is the average of a point dipole density located at Tc, and 
similarly for higher orders. The dipole must be understood as the average of 
the quantum mechanical dipole operator, as an exact analogy with the electric 
dipole of classical electrodynamics cannot be made. The density can thus be 
viewed as a collection of point multipoles, located at the ccntroid of each wave 
packet. The microscopic distribution of spin is important at the molecular level, 
but at the macroscopic level the effect of this molecular distribution is replaced 
by a sum of multipoles. Since the center of spin is different from the center of 
charge, in principle all multipoles are present. 

Following a similar manipulation and using the Boltzmann equation, the 
torque density is re-expressed as: 

T'^(R,i) = p^'^(R,t)- V-P^'^(R,i) (13) 

with the torque monopole density 

p^^{R,t) = J J d^kd''rJ{r„-k,t){T^)6{R-r,) = J d^kf{T^)\r,=n, (14) 

and the torque dipole density 

W-{R,t)= J J d^kd^rcf{rc,Kt){{R-r,)T^)\r^=n = J ti3W|r,=R. (15) 

In analogy with the spin dipole, the torque dipole has been defined as p^'^ = 
((f — rc)T'^} |r^=R. The torque density is therefore also a sum of multipole 
moments, that is, the moments of a point spin source located at Vc- Even in the 
case when the center of (s'^} coincides with the center of charge, (f'^) may not be 
centered at Tc, with the result that the higher order terms in the torque density 
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are in general present. The second and higher terms of cancel exactly the 
analogous terms in the continuity equation which come from the current. 

Since only the gradient of the spin current appears in the equation of conti- 
nuity, in the expansion of the sampling operator we keep the leading term 

J'^(R,t) = J d'kf{vsnL=n- (16) 

Keeping terms to first order in (f — Tc), the current can be decomposed into the 
following: 

The convective term C^°" represents the spin being transported along with the 
wave packet 

C^-(R,t) = J J d^kd^rJ{r,,-k,t)r,{s-)S{-R-r,) = J d^kfc'''\,^=R, (18) 

while D''"' comes from the rate of change of the spin dipole, which has already 
been introduced. It has the form: 

D-(R,<) = y'dV^k=R (19) 

P'^°' is the torque dipole introduced above. The corresponding monopole term 
appears in the source term of the continuity equation, as will emerge below. 
The presence of the torque dipole here is to be contrasted with the absence of 
an analogous term in classical electrodynamics. There, an electric dipole arises 
from the placement of two charges a small distance from each other, but the 
charge itself is conserved. 

Finally, we come to the source term in ([5]). The first part, composed of the 
torque density, has already been discussed. The second term, denoted by , 
becomes, in the relaxation time approximation 

= = - / d'Hfo - mn, (20) 

where Tp is the momentum relaxation time and /o the equilibrium distribution, 
which is usually the Fernii-Dirac distribution function. 

Based on the continuity equation alone, there is some flexibility in defining 
the current and the source. In systems in which spin is conserved, the torque 
density becomes, to first order in (f — Tc), a pure divergence, which can be 
incorporated into a redefinition of the spin current. This current, henceforth 
referred to as the spin transport current, is only due to the convective and spin 
dipole contributions: 

3*" (R, t) = (R, t) + B'" (R, t) (21) 



9 



With respect to this spin transport current, the continuity equation takes the 
following form: 



In the steady state under a constant electric field, the distribution function is 
composed of an equilibrium part, independent of the field, and a non-equilibrium 
part, which is first order in the field. Henceforth, terms in the spin current and 
source which depend on the equilibrium distribution function will be referred 
to as intrinsic, whereas the terms depending on the nonequilibrium shift in the 
distribution will be referred to as extrinsic. For example, the integrand in Eq. 
(|16p can be decomposed into a zero order spin- velocity, v(s'^), where v is the 
usual group velocity of the band, and a first order correction. Therefore, there 
will be a contribution to the current from the non-equilibrium part of the dis- 
tribution and the zero order spin-velocity, which has been discussed extensively 
in previous work [29l [30j [SU [32j [33] . There will also be a contribution from 
the equilibrium distribution and the first order correction to the spin-velocity, 
which is referred to as the intrinsic contribution. In the wave packet formalism 
presented here this effect arises from the change in wave functions induced by 
the electric field, rather than from the change in distribution functions that is 
responsible for most conventional transport effects. The intrinsic spin current 
is calculated from (jl6p using the equilibrium distribution and the expectation 
values of the spin and spin dipole operators in a Bloch state perturbed to first 
order in E. 

In its turn, the source in ^ can be decomposed into intrinsic and extrinsic 
contributions. The present entry considers homogeneous systems, so that all the 
gradient terms vanish, and the torque density is simply /(f °'). The zeroth order 
contribution to this term is null, as the Bloch wave functions are eigenstates 
of the Hamiltonian. Thus, to first order in the electric field, we find that (f'^) 
is simply given by (eE/h) ■ (9(s'^)/9k). One is thus justified in replacing / 
by its equilibrium value /o, in which case this term is purely intrinsic. The 
second term in the source, J-"^ , which depends on the nonequilibrium shift in 
the distribution function, is entirely extrinsic. 

The extrinsic source term takes into account the effect of scattering, and is 
a term which usually appears in the equation of continuity. The intrinsic source 
accounts for the effect of all spin-nonconserving terms, and must be present even 
in a clean system, if the Hamiltonian contains spin-dependent contributions. In 
general, in addition to the rate of change of spin arising from the spin-dependent 
terms in the Hamiltonian, scattering processes may alter the orientation of the 
spin, with the result that any one spin component is not conserved, and the 
orientation of spins is randomized over a longer time period. For a uniform 
steady-state system, the current is constant and the intrinsic source term must 
vanish. However, near the boundary of the system, or at an interface with a 
different semiconductor with (for example) weaker spin-orbit interactions, the 
spin current driven by an electric field will vary spatially and T must reach a 
non-zero value. 




(22) 
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Let us take a closer look at the spin dipole and torque dipole, which are 
seen to be the main mechanisms responsible for generating the spin current. 
Because of its narrow distribution in k, the mean spin of the wave packet is 
{wi\s'^\wi) = {ui\s'^\ui), where it is understood that the wave vector of the Bloch 
function is set at kc and s'^ is an arbitrary projection of the spin vector operator. 
The spin dipole of the wave packet, defined relative to the charge center of the 
wave packet is given, in terms of Bloch functions, by the expression: 

pr = ^[("«l^n^> - (^1^1"^)] - {n^\^^){u.\s^\u..) (23) 

Interestingly, the spin dipole is independent of the wave packet width. The 
expression is also invariant under a local gauge transformation, in the sense 
that if \ui) is modified by a phase factor e*"'-'''' the spin dipole is unchanged. 

The torque dipole term has a special interpretation in the case of spin trans- 
port. The rate of change of spin is equivalent to a torque, and the torque dipole 
represents the moment exerted by this torque about the center of the wave 
packet. The semiclassical expression for the torque moment is 

pr = ^[(^^i^n^) - (^i^'^i"^)] - {^^\^^)Msnu,). (24) 

The torque moment has the same gauge invariance properties as the spin dipole, 
and like the spin dipole it also does not depend on the wave packet width. 
It is important to note that the spin current, JT*^, can be simplified to: 

J%^,t)= j d'kfXv{u,\rw\u,), (25) 

which is the semiclassical equivalent of the Kubo formula for spin currents. 

3.2 Density matrix picture of spin transport 

Semiclassical theory provides a straightforward, intuitive picture of the way 
spin currents arise in the course of carrier dynamics in an electric field. The 
theory was developed for independent bands. It turns out that interband co- 
herence arising from scattering is crucial in spin transport, and is difficult to 
treat semiclassically. Although the semiclassical theory can be generalized to 
multiple bands [491 150j , it is more instructive to examine spin transport from a 
different point of view that is closer in outlook to the philosophy underlying the 
Kubo formula (with which the semiclassical theory agrees.) This will shed some 
light on additional issues, such as the relationships between spin currents and 
spin precession, between spin currents and spin densities, the complex effect of 
disorder and the vanishing of spin current in certain systems. 

A large, uniform system of non-interacting spin- 1/2 electrons is represented 
by a one-particle density operator p. The expectation value of an observable 
represented by a Hermitian operator O is given by tr(/50) and p satisfies the 
quantum Liouville equation 

| + l[H + a,^1^0. (26) 
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The Liouville equation is projected onto a set of time-independent states of 
definite wave vector {|fes)}, which are not assumed to be eigenstates of the 
Hamiltonian H . The matrix elements of p in this basis will be written as pkk' = 
Pkk' = {ks\p\k's'). Spin indices are not shown exphcitly, and pkk' is a matrix 
in spin space, referred to as the density matrix. In this work we require the 
expectation values of operators which are diagonal in wave vector, and will thus 
require the part of the density matrix diagonal in wave vector, pkk = fk = nk^ + 
Sk- In the presence of a constant uniform electric field fk = fok + fsk, where 
the equilibrium density matrix /ofc is given by the Fermi-Dirac distribution, 
and the correction fsk is due to the E. We subdivide fok = n^k 1 + S^k and 
fsk = T^Bfc 1 + SEk- The spin-dependent part of the nonequilibrium correction 
to the density matrix Ssk is interpreted as the spin density induced by E. The 
equations governing the time evolution of UEk and Ssk is |51| 



duEk ~, I s e.E dnok 

H - [-tik, JEk\ + Jo (JEk) — • —gj^ Js [nEk) = ^Ek, 



(27) 



where the scalar part of the scattering operator Jq and its spin-dependent part 
Js have been defined in Ref. |51| . The equation for nEk has the well-known 
solution UEk = (eETp/h) ■ (driok/dk), in other words, UEk describes the shift 
of the Fermi sphere in the presence of the electric field E, with the momentum 
relaxation time Tp. It is seen from Eq. (P7|) that spin-dependent scattering gives 
rise to a renormalization of the driving term in the equation for SEk- This 
renormalization has no analog in charge transport. 

We need to find the expectation value of the spin current operator defined 
in the introduction. In the systems under study the spin current operator can be 
written as = hhs" /m* + {l/^K) dQ'^/dht. We need to determine Ssk- To 
this end we remember that an electron spin at wave vector k precesses about an 
effective magnetic field flk- The spin can be resolved into components parallel 
and perpendicular to fife. In the course of spin precession the component of 
the spin parallel to flk is conserved, while the perpendicular component is 
continually changing. Corresponding to this decomposition of the spin is an 
analogous decomposition of the spin distribution S Ek into a part representing 
conserved spin and a part representing precessing spin, denoted by SEk\\ and 
SEk± respectively. There is an analogous decomposition of the source on the 
RHS of Eq. P7)) into T,Ek\\ and ^Ek±- This decomposition is carried out by 
introducing projection operators P|| and Pi as described in Ref. |51| . giving for 
SEk\\ and SEk± in the weak momentum scattering limit 



dS 



Ek\\ 



+ P\lJoiSEk) = ^Ek\u (28a) 



dt 

+ ^ [Hk,SEk±] = ^Ek± - P±Jo {SEk). (28b) 
Equation (I28b|) shows that scattering mixes the distributions of conserved and 
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precessing spins. This is so because when one spin at wave vector k and pre- 
cessing about flk is scattered to wave vector k' and precesses about Clk', its 
conserved component changes, a process which alters the distributions of con- 
served and precessing spin. Equations (pS)) can be solved straightforwardly if 
one assumes the impurity potential to be short-ranged, obtaining |51| 

SskW = ^Ek\\ Tp + P|| (1 - P\\)~^^Ek\\ Tp, (29a) 

_ flk X i'^Ek±Tp + P± SEk\\) ■ CTTp {T,Ek±Tp + P± SEk\\\ ,s 

""^^ " 2h{i + niT^/h^) 1 + niT^/h^ ^^^^^ 

The correction SEk\\ does not give rise to a spin current. Inspection of Eq. (|29ap 
shows that integrals of the form J dO Sf[ SEk\\ contain an odd number of pow- 
ers of k and are therefore zero. It can, however, give rise to a nonequilibrium 
spin density, since integrals of the form J d9 s" SEk\\ contain an even number of 
powers of k and may be nonzero. Similarly SEk± does not lead to a nonequi- 
librium spin density. The expectation value of the spin operator yields integrals 
of the form J dO s'^ 5'^^^ , which involve odd numbers of powers of k and are 
therefore zero. This term does, however, give rise to nonzero spin currents, 
since integrals if the form J dO SEk± contain an even numbers of powers 
of k and may be nonzero. Therefore, in the absence of spin-orbit coupling in 
the scattering potential, nonequilibrium spin currents arise from spin precession 
(as outlined by Sinova et al. [35 ), and nonequilibrium spin densities from the 
absence of spin precession. The dominant contribution to the nonequilibrium 
spin density in an electric field exists because in the course of spin precession 
a component of each individual spin is preserved. For an electron with wave 
vector fc, this spin component is parallel to fife. In equilibrium the average 
of these conserved components is zero. When an electric field is applied the 
Fermi surface is shifted and the average of the conserved spin components may 
be nonzero, as illustrated in Fig. 1. This argument explains why the nonequi- 
librium spin density (x t^^ and requires scattering to balance the drift of the 
Fermi surface. Although spin densities in electric fields require band structure 
spin-orbit interactions and therefore spin precession, the dominant contribution 
arises as a result of the absence of spin precession. 

Systems in which is linear in k are special, in that the spin current as 
defined in this section vanishes [39l|40l|4Tl|Mlllllllll3[46]. This is because 
of the renormalization of the driving term on the RHS of Eq. (I28b|) for SEk±, 
in other words because of scattering from the conserved spin distribution to the 
precessing spin distribution. In Eq. (j29bp it is also clear that if T^EkJ^Tp + 
P-L SEk\\ vanishes, then all the contributions to SEk± also vanish. Since SEk\\ 
effectively represents a steady-state spin density, we see that the presence of 
this spin density tends to diminish the spin current. In systems with energy 
dispersion linear in k it cancels the spin current completely. 
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4 Future Directions 

Whereas the community appears to be in agreement that spin currents exist and 
are measurable, many questions remain unanswered. Theoreticahy, intrinsic and 
extrinsic effects (such as due to skew scattering and side jump) have not been 
studied on the same footing for an arbitrary form of band structure spin-orbit 
interactions. The relative magnitude of intrinsic and extrinsic spin currents in 
such a general system remains to be determined. Also, different definitions of 
the spin current give results that often differ by a sign [H [5] . The relationship 
between spin current and spin accumulation at the boundary is not clear, again 
thanks to the non-conservation of spin. It appears that what happens at the 
boundary is sensitive to the type of boundary conditions assumed. Thus so 
far as quantitative interpretation of experimental data is concerned, theory has 
some way to go. 

Despite tremendous progress, experiment is still searching for a reliable way 
to measure, as opposed to detect, spin currents directly. Practically, the ques- 
tion of what to do with spin once it has been transported/generated remains. 
The revolutionary electronic device that harnesses spin currents for a practical 
purpose remains to be made, and the challenge of its design confronts experi- 
mentalists and theorists alike. 

The research at Argonne National Laboratory was supported by the US 
Department of Energy, Office of Science, Office of Basic Energy Sciences, under 
Contract No. DE-AC02-06CH11357. 
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